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Abstract 

Both the Sikorsky Finite-Element Airframe 
Vibration Analysis Program (FRAN /Vibration Analy- 
sis) and the NASA Structural Analysis Program 
(NASTRAN) have been correlated vith data taken 
in full-scale vibration tests of a modified CH-53A 
helicopter. With these programs the freq.uencies 
of fundamental fuselage bending and transmission 
modes can be predicted to an average accuracy of 
three percent with corresponding accuracy in 
system mode shapes . 

The correlation achieved with each program 
provides the material for a discussion of modeling 
techniques developed for general application to 
finite-element dynamic analyses of helicopter 
airframes. Included are the selection of static 
and dynamic degrees of freedom, cockpit structural 
modeling, and the extent of flexible- frame model- 
ing in the transmission support region and in the 
vicinity of large cut-outs. The sensitivity of 
predicted resiilts to these modeling assixmpticaas 
is discussed. 

Introduction 

Helicopter vibration and resulting aircraft 
vibratory stress can lead to costly schedule 
slippages as well as to problems in field service 
maintenance and aircraft availability . At the 
core of vibration control technology is the require- 
ment to design the helicopter structure to minimize 
structural response to rotor excitations. Both the 
complexity of the structure and the increasingly 
stringent mission and vibration control specifica- 
tions demand development of airframe structural 
vibration analyses that can be used rapidly and 
economically to evaluate and eliminate vibration 
problems during the prelimineiry design phase of 
helicopters. 

The complex helicopter structure consists of 
sections that differ considerably in structural 
arrangement and load carrying requirements. These 
sections include the cockpit, cabin, tail cone, 
and tail rotor pylon. In addition, large fuselage 
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cut-outs and concentrated masses such as the trans- 
mission, main rotor, and tail rotor, which are 
unique to helicopters, play a major role in con- 
trolling vibrations. 

Although advanced analytical methods based on 
finite-element techniques have been developed for 
studying the vibration characteristics of con 5 >lex 
structures, a detailed correlation of such methods 
with test data is not available in the general 
literature. Further, little information is avail- 
able on the accuracy of various modeling asstimp- 
tions that might be made to reduce the cost and 
time of applying these vibration analyses. 

As a resTilt a research project was establisb- 
by Naval Air Systems Command with Sikorsky Aircraft 
to: 

a) Determine the accuracy of the Sikorsky 
Finite-Element Airframe Vibration 
Analysis in predicting the vibration 
characteristics of complex helicopter 
airframe structures. 

and 

b) Develop and evaluate general helicopter 
dynamic modeling techniques that could 
be used to provide accurate estimates 

of vehicle dynamic characteristics vhxle 
at the same time minimizing the com- 
plexity and cost of the analysis. 

Due to the increased usage of NASTRAN 
throughout the industry as well as the efficiency 
resulting from employing a single analytical sys- 
tem for both static and dynamic analyses, a par- 
allel correlation study using NASTRAN has been 
performed. The results of these correlation 
studies are the subject of this paper. 

Phase I - Stripped Vehicle 

Test Vehicle 

At the initiation of this effort, the phi- 
losophy guiding the development of modeling tech- 
niques was based upon the concept of gradually in- 
creasing the complexity of the analytical repre- 
sentation. It was decided that the first 
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correlation study would be conducted on an air- 
craft stripped of all appendages. It was believed 
that the modeling techniques for representing air- 
frame response characteristics could be identified 
and developed most easily in this manner. Then, 
as various appendages were added to the basic 
vehicle, only the modeling techniques required for 
the structure or masses added need be developed. 

The vehicle used in this test and correla- 
tion study was the CH-53A Tie Down Aircraft, Ve- 
hicle designation number 6l3. A general arrange- 
ment of the structure is illustrated in Fig\ire 1. 
For initial correlation, all appendages were re- 
moved. These included the nose gear, main landing 
gear, main landing gear sponsons, fuel sponsons, 
tail pylon aft of the fold hinge, tail rotor and 
associated gear boxes, engines, cargo ramp door, 
horizontal stabilizer, and all remaining electri- 
cal and hydraulic systems. The main rotor shaft 
and all gears were removed from the main transmis- 
sion housing and only the housing itself was re- 
tained for the test configuration. 

Testing 

The ground test facility employed to estab- 
lish the dynamic characteristics of the test vehi- 
cle was a bungee suspension system that simulates 
a free-free condition, a rotorhead-mounted uni- 
directional shaker, and the Sikorsky shake test 
instrumentation console. Instrumentation con- 
sisted of ik fixed and 10 roving accelerometers. 

A complete description of the test apparatus and 
the instrumentation is provided in Reference 1. 

All accelerometer signals and the reference 
shaker contactor signal were transmitted to the 
console. The signals were processed automatically 
by the console resulting in a calculation of the 
in-phase and quadrature components of the acceler- 
ations. The accelerations were then normalized to 
the magnitude of the shaker force at the particu- 
lar frequency. As frequency was varied, the re- 
sulting response of each accelerometer was record- 
ed on a XYY* plotter, Figure 2, as g* s/1000 lbs. 
versus frequency. 

Ideally, a fuselage mode can be identified 
by a peak in the quadrature response and a simul- 
taneous zero crossing of the in-phase response. 
Once a mode is located, all quadrature responses 
at this frequency can be recorded to define the 
mode shape. The modes defined in this manner from 
the shake tests are listed in the left-hand column 
of Table I. It should be noted that this tech- 
nique is applied more easily at lower frequencies, 
where sufficient modal separation exists so that 
the forced response in the vicinity of a resonance 
is dominated by a single mode. As shown in Figure 
2, the mode shapes at higher frequencies must be 
extracted from the coupled response of many modes. 

Analysis and Correlation 

The shake test data indicated that the 
natural modes of vibration of a helicopter can be 
categorized as beam-like modes controlled by 


overall fuselage characteristics (e.g., length, 
depth, overall bending stiffness, mass distribu- 
tion, etc.) and those controlled by the transmis- 
sion support structure. Therefore, the overall 
helicopter structure was modeled utilizing three 
modules : 

1) center section including the transmission 

support region 

2) forward fuselage and cockpit 
and 

3) aft fuselage and tail. 

The center section was modeled in greatest detail 
by applying finite-element techniques. The struc- 
tural characteristics of the forward and aft fuse- 
lage were derived from beam theory. These equiva- 
lent beams were located at the neuLral axis of the 
airframe section and were assigned the bending and 
torsional properties of the total section. The 
beam models of the forward and aft fuselage were 
cantilevered from rigid frames at the respective 
forward and aft ends of the center section. Figure 
3. The influence coefficients of these beams with 
respect to their cantilevered ends were then com- 
bined with the influence coefficient matrix of the 
remaining structure. 

The Phase I correlation was performed using 
the Sikorsky Finite -Element Airframe Vibration 
Analysis ( FRAN /Vibr at ion Analysis) . This analysis 
consists of two programs: PPFRAN and a 200 dynamic- 

degree-of-freedom eigenvalue/eigenvector extraction 
procedure. PPFRAN is derived from the IBM/MIT 
Frame Structural Analysis Program, FRAN (Reference 
2), a stiffness method, finite-element analysis 
limited to two types of elements, namely bending 
elements (bars) and axial elements (rods). This 
limitation necessitated further development of FRAN 
for application to stressed skin structures. This 
development consists of the addition of pre- and 
post- operative procedures linked to FRAN. In the 
pre-operative procedure (Pre-FRAN), the fuselage 
skin is transformed into equivalent rod elements. 
This transformation is developed by satisfying the 
criterion that the internal energy of the skin 
structure under an arbitrary set of loads be the 
same as that of the transformed structure under the 
same set of loads. The post-operative procedure 
(Post-FRAN) extracts the influence coefficient 
matrix corresponding to the selected dynamic 
degrees of freedom. A detailed description of the 
FRAN /Vibration Analysis is provided in Reference 1. 

The elements used to represent the airframe 
structure are: 

1) bending (bar) elements for fuselage frames 

and for the nose and tail beams 

2) axial (rod) elements for the stringers 
and 

3) equivalent, diagonal rod elements for skin 

panels . 
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For dynamic analysis, the structure is as- 
sumed to be unbuckled, so that all skin panels are 
considered fully effective in resisting axial ^ 
loads* Thus, the total axial area of each skin 
panel is lumped with the areas of adjacent string- 
ers. 


During Phase I correlation, three modeling 
parameters were varied: the number of bays over 

which the finite-element (flexible -frame) model 
extends (Figure U), the number of nodes per frame 
(niamber of stringers), and the number of dynamic 
degrees of freedcan assigned to each frame (Figure 
5). The results of the correlation are presented 
in Table I, which shows the sensitivity of the 
analysis to each of the above parameters and the 
accuracy of the predicted frequencies and mode 
shapes. The criteria for establishing the level 
of mode shape correlation are: 

E (Excellent) - Correct number of nodes, nodes 

less th8in 2.5 percent of fuselage 
length from measirred location, 
local modal amplitudes within 20 
percent of test values. 

G (Good) - Correct number of nodes, nodes 

less than 2.5 percent of fuselage 
length from measured location, 
difference between actual and pre- 
dicted local modal amplitudes ex- 
ceeds ±20 percent of test values. 

F (Fair) - Correct number of nodes, nodes 

more than 2.5 percent of fuselage 
length from measured location, 
difference between actual and pre- 
dicted local modal amplitudes ex- 
ceeds ±20 percent of test values, 

P (Poor) - Incorrect number of nodes, nodes 

located improperly, difference be- 
tween actual and predicted local 
modal amplitudes exceeds ±20 per- 
cent of test values. 

A comparison of the 30- and 60-strlnger anal- 
yses indicates that there is no change in the re- 
sults when modeling the structure with half the 
number of actual stringers. In addition a compari- 
son of results obtained with the basic and reduced 
dynamic degree of freedom allocation indicates that 
no more than l6 dynamic degrees of freedom per 
frame are required for dynamic modeling. 

Although mode shape correlation res\ilting 
from the analysis in the frequency range of inter- 
est (below 1500 cpm) is encouraging, see Table I, 
the absence of a representative mass distribution 
made the analysis overly sensitive to certain mod- 
eling assumptions. This sensitivity appears to 
accoxmt for the less than satisfactory frequency 
correlation. For example, the frequency of the 
transmission pitch mode is normally controlled by 
the mass of the fully assembled transmission and 
the properties of the structure in the transmission 
support region. In the absence of a mass distribu- 
tion representative of a fully assembled vehicle, 
however, any element of the structure and any 


lunsped mass can contribute significantly to the 
control of the dynamic characteristics. In this 
case, the analytical representation appears to be 
too stiff because of the beam model used for the 
fuselage forward of F.S. 262, which constrains the 
upper and lower decks to deform equally. Thxs 
constraint is not imposed by the actual structure. 

A ccaaparison of the results of the 9- and l8-bay 
analyses indicates that due to the local nature of 
the transmission pitch mode, extension of the 
flexible-frame model aft beyond the limit of the 
9-bay model has no significant effect on the pre- 
diction of this mode. 

The poor frequency correlation for the first 
lateral bending mode persisted throughout this 
phase of correlation. This mode was characterized 
by differential shearing of the upper and lower 
decks of the aft cabin, Fig\ire 7- The 6-bay and 
9-bay flexible- frame model represented most of this 
structure experienceing the differential shearing 
as a beam capable of only bending and torsion. 

This overly constrained model resiilted in predicted 
frequencies substantially higher than test values. 
Extending the flexible-frame representation to 18- 
bays appears to be the solution. However, size 
limitations in PFFRAN required that the l8— hay 
flexible-frame model be generated in two 9-bay sub- 
structures, married at a rigid intermediate frame 
at F.S. Figure 3. Although the extended model 

improved the correlation of the first lateral bend- 
ing mode, absence of a representative mass distribu- 
tion again appears to make the model overly sensi- 
tive to the presence of the rigid frame at F.S. ^*^^2. 
This accounted for the remaining difference between 
test and analysis. 

Many of the higher frequency modes are con- 
trolled by the structure in the area of the rear 
cargo ramp. This accounts for the failure to pre- 
dict the Transmission Vertical mode until the flex- 
ible-frame model was extended into the ramp area, 
see Table I. Although this extension of the model 
improved correlation, the high frequency modes 
above 1500 cpm are difficult to identify analytical- 
ly due to the coupling of overall fuselage modes 
with local frame modes. This difficulty is com- 
pounded in this investigation, because the fre- 
quencies of the basic fuselage modes are raised due 
to the stripped condition of the vehicle, while 
frequencies of the local frame modes are lowered 
due to the lumped-mass modeling used to represent 
each frame. Tests of a more representatively load- 
ed fuselage can be expected to minimize the problem 
of mode identification. 

From the results of this phase of the corre- 
lation, it is concluded: 

1) The selection of static degrees of freedom in 
the flexible frame model can be based on a 
structural model that contains stringers num- 
bering one half the number of actual 
stringers . 

2 ) 
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No more than sixteen dynamic degrees of free- 
dom on each flexible frame are required for 
dynamic analysis. The typical location of 



these degrees of freedom is illustrated in 
Figure 5. 

3) Transmission modes can he predicted by a 

flexible- frame representation of the trans- 
mission support region extending about I.5 
transmission lengths forward and aft of the 
corresponding transmission supports, about 
9 bays. If the vehicle contains large cut- 
outs, such as the cargo ramp of the test 
vehicle, the flexible-frame model should ex- 
tend through this region as well. 

PHASE II - BALLASTED VEHICLE 

Testing 

Shake tests were performed after adding bal- 
last to provide a more realistic representation of 
a helicopter mass distribution. Figure 6. At the 
transmission mountiag plate, two lead blocks hav- 
ing a total weight of ^570 pounds were mounted so 
that the mass and pitching moment of inertia of 
the simulated transmission and rotor head approx- 
imated that of the actual CH-53A. At the tail, a 
1500-pound block was mounted to simulate the re- 
moved tail pylon, stabilizer, and tail rotor. At 
the nose, a 3000-pound block was mounted on the 
nose gear trunnion fitting to balance the vehicle. 

The natural modes of vibration identified by 
shake tests are listed in Table II along with the 
frequencies measured during Phase I. Not only did 
the ballast succeed in lowering the fuselage modes 
into a frequency range more representative of that 
encountered on a fully assembled aircraft, but ad- 
ditional modes were also identified that etre 
strongly controlled by the ballast. In fact, 
these modes were identified as local modes of the 
ballast blocks themselves. Due to the complex 
structural nature of the ballast. Figure 6, these 
appendages did not lend themselves to simple ana- 
lytical representations. Therefore, the flexi- 
bility of each ballast block was measured by in- 
strumenting both the block and the adjacent air- 
frame structure and then measuring the accelera- 
tions occurring at both locations near the modal 
frequencies of interest. The mass of each ballast 
block and its absolute acceleration resulted in a 
force which produced the relative motion between 
the two instrumented parts. The empirically de- 
fined flexibilities of the ballast were then used 
in the dynamic model. 

Analysis and Correlation 

The modeling techniques developed in Phase I 
of this study were applied to both the FRAN/Vibra- 
tion Analysis and NASTRAN. 

The finite-element model analyzed in Phase 
II was identical to the l8-bay model analyzed in 
Phase I, except for adding the mass and structural 
characteristics of the ballast blocks. The FRAN 
model was formed with rod and bar elements, as 
discussed previously, while the NASTRAN model used 
CROD, CBAR, and CSHEAR elements (Reference 3). As 
before, all skin panels were asstuned fully 


effective in reacting axial load and this effective 
area was lumped into the adjacent stringers. 

Including ballast to replace removed appenda- 
ges resulted in a substantial improvement in the 
correlation, particularly in frequency prediction 
as shown in Table III, Significantly, ballast 
eliminated the difficulties identified as sensitivi- 
ty to modeling assumptions and local frame modes in 
the absence of representative mass distributions. 

The average error in predicting the frequencies of 
fundamental fuselage bending modes and the trans- 
mission pitch mode was 3-W for both the FRAN /Vibra- 
tion Analysis and NASTRAN. In addition the shape 
correlation for these modes varied from good to ex- 
cellent. The analyses also were able to predict 
accurately the significant changes in the charac- 
teristics of the fuselage and transmission modes 
resulting from the addition of the ballast. Figures 
7, 8 and 9« To achieve this degree of correlation, 
modeling of the ballast flexibilities was required. 
This modeling was successfully accomplished in the 
vertical/pitch direction. Figure 10, but did not 
prove successful in the lateral/torsion direction, 
Figure 11. The contrast between these two results 
establishes the ability of finite -element analyses 
to predict accurately the characteristics of fuse- 
lage and transmission modes when the structxiral 
data base is defined with sufficient accuracy. Fur- 
ther improvement in the correlation could have been 
achieved if a more detailed definition of the bal- 
last flexibilities had been acquired from measure- 
ments of static deflections. 

Reasonable success has been achieved in pre- 
dicting higher frequency, ramp-controlled modes. 
Figures 12 and 13. However, some margin does exist 
for further improvements in shape and frequency 
prediction. From the steuidpoint of modeling, it 
appears that the 200 dynamic degree of freedom 
limit established in this study is inadequate for 
predicting the shell-type modes of the cargo ramp 
structure. In addition, the test procedure em- 
ployed, namely the use of a single rotorhead 
shaker, does not provide a means of uncoupling the 
forced response characteristics of modes at the 
higher frequencies. Figure 2. 

Conclusions 

1. Finite element analyses can predict accurate- 
ly the frequencies and mode shapes of complex 
helicopter structiires, provided the structur- 
al data base is defined accurately, 

2. Complete stripping of a vehicle for correla- 
tion purposes may make the analysis overly 
sensitive to normally minor modeling assump- 
tions. 

3. Significant changes can be predicted accurate- 
ly in the character and frequency of fuselage 
and transmission modes due to changes in mass 
distributions and structural characteristics. 

The modeling techniques established by this 
study can be used during aircraft design re- 
gardless of the finite-element analytical 
system being used. 
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Figure 2 Typical Vertical Response to Vertical Excitation 


TABLE I 

SHAKE TEST CORRELATION SUMMARY 



l8 Bay 
30 Stringer 
Reduced DOF 

9 Bay 

30 Stringer 
Reduced DOF 

6 Bay 

30 Stringer 
Reduced DOF 

30 Stringer and 
60 Stringer 
Basic DOF 

■ 

Mode 

Freq. 

(CPM) 

Freq. 

Error 

Shape 

Freq. 

Error 

Shape 

Freq. 

Error 

Shape 

Freq. 

Error 

H 

1st Lateral 

910 

1207 

33% 

G 

IU66 

60^ 

P 

li+35 

58% 

P 

1I4UO 

5Q% 

p 

1st Vertical 

1155 

1175 

2% 

E 

1282 

11^ 

E 

I2U2 

8 % 

E 

12l4l 

8 % 

E 

XSSN Pitch 

1^90 

1709 

13% 

E 

1710 

13^ 

E 

171+8 

n% 

E 

1758 

11 % 

E 

2nd Vertical 

1950 

2150 

10% 

G 

2390 

22% 

F 

2505 

28 % 

F 

2577 

32 % 

F 

XSSN Roll 

2000 

2U05 

20% 

P 

2870 

U3% 

P 

2900 

h^% 

P 

289I+ 

1+5^ 

P 

XSSN Vertical 

2150 ' 

2250 

k% 

F 










Torsion 

2300 

2763 

20% 

F 

2I428 

6 % 

P/F 

2k22 

6 % 

P/F 

2I+I45 

(y% 

P/F 
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Figure 6b Tail Ballast 
Installation 


Figure 6c Nose Ballast 
Installation 


Figure 6 Phase II Ballast Installations 
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TABLE II - SHAKE TEST FREQUENCIES 


MODE 

PHASE II 

TEST FREQUENCY CM 

1st Vertical Bending 

Uo 

1st Lateral Bending 

615 

Transmission Pitch 

IkQ 

Nose Block Lateral /Roll 

930 

Nose Block Vertical/Coupled 

970 

Forward Cabin/Nose Block Lateral 

990 

Nose Block Vertical 

1050 

Second Vertical Bending 

1290 

Torsion 

1310 

Tremsmission/Ramp Vertical Bending 

ll*25 

Ramp Vertical Bending 

l6l+0 


PHASE I 
TEST FREQMfCY 

1155 

910 

lJ*90 


1950 

2300 

2150 


MODE 

1st Vertical Bending 

Transmission Pitch 

Nose Block Vertical/ 
Transmission Pitch 

Nose Block Vertical 

Second Vertical 

Transmission Vertical/ 
Ramp Vertical 

Ramp Vertical 


1st Lateral Bending 
Nose Block Lateral/Roll 

Forward Cabin Lateral/ 
Nose Block Lateral 

Torsion 


TABLE III 

PHASE II CORRELATION SUMMARY 
VERTICAL/PITCH MODES 


Frequency 


Test 

FRAN 

Error 

Shape 

1*1»0 

k3S 

0? 

E 

7lt0 

751 

1.5? 

G 

970 

933 

w 

G 

1050 

10U3 

1? 

F 

1290 

1523 

H 

00 

F 

lk25 

1563 

10? 

F/G 

l6U0 

139>* 

15? 

P/F 


T.A'miAL/TORSION MODES 

615 

659 

T? 

G 

930 

735 

21% 

P 

990 

858 

13? 

P 

1310 

1601 

22? 

P 


NASTRAB 

Error 

Shape ! 

U53 

3 ? 

E 

i 

785 

6 ? 

G 

956 

1 . 5 ? 

G 

1063 

1 ? 

F j 

1608 

25? 

F ' 

I8U3 

29? 

F/G 

1355 

17 ? 

P/F 

1 

595 

3 ? 

1 

1 

i 

G 

812 

13 ? 

P 

970 

2 ? 

P 

1325 

1 ? 

P 
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